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ABSTRACT 

A  method  is  described  for  the  introduction  into  the  Fortran  63 
language  of  a  special  operator  for  differentiation  which  would  be 
executed  when  encountered  during  compilation.     The  original 
expression  would  be  extracted  from  the  Fortran  equation  and 
replaced  by  its  symbolic  derivative;  the  translation  process  would 
then  proceed  in  the  normal  manner. 

A  model  algebraic  translator,  written  in  Fortran  63,  is 
described  and  test  results  are  presented. 
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I.       INTRODUCTION 

Hanson    and  others  before  him  recognized  the  need  for  a 
compiler  routine  which  would  accept  an  algebraic  expression  and 
print-out  a  symbolic  representation  of  the  derivative  of  that 
expression  in  addition  to  generating  the  necessary  machine 
symbolic  macros  required  for  subsequent  numerical  evaluation  by 
the  computer.    This  thesis  presents  a  method  to  accomplish  those 
goals  in  the  Fortran  63  language.  \_3J 

The  method  described  herein  assumes  the 'utilization  of 
reverse  Polish  notation2  and  a  left  to  right  scan  during  the 
compilation  process. 

"DX"  is  defined  as  a  special  operator,  signifying 
differentiation.    This  operator  is  somewhat  analogous  in  design 
and  use  to  the  normal  Fortran  "Function",  e.g.  ,  SINF,  which 
appears  to  be  treated  both  as  an  operator  and  an  operand  during 
its  life  cycle.    Certainly,  SINF  signifies  and  initiates  an 
operation.    Yet,  during  the  translation  process,  the  function  is 
treated  as  an  operand  and  is  paired  with  a  special  operator, 


Hanson,  J.  W.  ,  J.  S.  Caviness  and  C.  Joseph.     "Analytic 
Differentiation  by  Computer",  Communications,  ACM,  Vol.  5, 
No.   6,  June  1962:  349. 

Hamblin,  C.  L.  ,   "Translation  to  and  from  Polish  notation", 
Computer  Tournal  5 ,  3,  October  1962:    210-213. 


usually  the  period  symbol.     This  similarity  in  the  utilization  of  the 
"DX"  will  be  evident.     However,   since  there  will  also  be 
dissimilarities,  the  authors  have  chosen  to  refer  to  the 
differentiation  operator,  DX,  as  a  "dual-operator"  for  clarity  of 
exposition. 

The  routine  for  translation  from  Fortran  notation  to  reverse 
Polish  notation  requires  the  recognition  of  the  dual-operator,  DX, 
and  a  special  subroutine  for  translation  of  the  dual-operator  and 
its  operand.    For  the  purpose  of  describing  this  procedure,  a 
simplified  routine  for  translation  to  reverse  Polish  is  presented. 

The  next  operation  of  translation  of  reverse  Polish  notation 
into  machine  symbolic  language,  likewise,  requires  recognition 
of  the  dual-operator,  DX,  and  a  special  subroutine  for  execution 
of  the  derivative  operator.    This  subroutine  is  designed  to 
extract  the  derivative  operator  and  its  operand  from  the  Polish 
string,  execute  the  differentiation,  place  the  symbolic  result 
into  the  Polish  string,  and  return  control  to  the  primary  routine 
at  the  original  point  of  entry. 
II.     HISTORICAL  BACKGROUND 

Hanson    presented  a  program  which  would  accomplish  the 
analytic  differentiation  of  a  mathematical  expression  on  a  digital 

3Hanson,  J.  W.  ,  et  al.  ,  op.  cit.  '/ 349-355. 


computer.     He  first  described  the  methods  and  accomplishments 
which  had  taken  place  in  this  field  prior  to  the  writing  of  his 
paper.     His  approach  was  to  take  an  input  mathematical  expression 
and  covert  it  into  triples.     He  then  took  the  differentiated 
expression  of  each  of  these  triples  incorporating  the  Ershov 
algorithm     in  the  differentiation. 

In  February  1964,   R.  E.  Wengert     wrote  a  paper  which 
developed  a  program  to  give  a  numerical  evaluation  of  a 
derivative  without  developing  the  analytic  expression.     His 
method  was  to  .set  up  a  library  of  elementary  function  subroutines. 
This  would  work  only  for  the  simple  routines ,  and  would  involve 

the  programmer  in  using  the  CALL  statements  in  his  program  to 

fi  7 

arrive  at  the  solution  to  his  problem.     Wilkens"  and  Bellman    also 

propose  the  use  of  library  subroutines  to  determine  the  value  of 

a  derivative. 


4Ibid. 

^Wenger,  R.  E.  ,   "A  Simple  Automatic  Derivative  Evaluation 
Program",  Communications,  ACM,  Vol.   7,  No.  8,  August  1964: 

463-464. 

i 

Wilkins,  R.  D.  ,   "Investigation  of  a  New  Analytical  Method 
for  Numerical  Derivative  Evaluation",  Communications ,  ACM, 
Vol.   7,  No.  8,  August  1964:    465-471. 

7 
Bellman,  R.  E.  ,  H.   Kagiwada,  R.  E.   Kalaba;  "Wengert's 

Numerical  Method  for  Partial  Derivatives,  Orbit  Determination 

and  Quasilinearization"  ,  Communications  ,ACM ,  Vol.  8,  No.  4, 

April  1965:  231-232. 
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Slagle    and  several  others  who  have  written  on  the  subject 
have  used  list  processing  to  accomplish  the  taking  of  the 
derivative.    Of  course,  this  process  cannot  be  used  by  any 
computer  which  does  not  have  the  list  processing  capability. 

However,  for  additional  comment  on  the  disadvantages  of  this 

9 
process,   see  Smith.       The  paper  by  Peter  J,  Smith  was  received 

too  late  for  a  true  evaluation  by  the  authors.     However,  he 

proposed  differentiating  in  symbolic  form  without  the  use  of  list 

processing,  or  the  requiring  of  the  programmer  to  write  subroutines 

as  a  part  of  the  program.    There  is,  however,  no  place  in  Smith's 

paper  any  mention  of  getting  a  numerical  answer  to  the  problem 

as  well  as  the  symbolic  representation  of  the  derivative. 

A  means  is  needed  whereby  the  symbolic  representation  of 

the  derivative  is  available,  as  well  as  the  numerical  value  for 

that  derivative. 


^Smith,  P.  J.  ,  Symbolic  Derivatives  without  List  Processing, 
Subroutines,  or  Recursion,  Ballistic  Research  Laboratories 
Memorandum  Report  No.   1630,  February  1965. 

9Ibid. 


III.    METHOD 

1 .    Definition  of  an  operator  for  differentiation . 

An  operator  usually  must  perform  two  functions:    identify  its 
operands  and  execute  its  designed  operation  upon  those  operands. 
In  addition,  a  derivative  operator  must  identify  the  variable (s) 
of  differentiation.    To  facilitate  the  performance  of  those  functions, 
the  authors  propose  the  adaptation  of  the  dual  nature  accorded 
the  "Function"  within  the  Fortran  language. 

For  use  in  the  Fortran  language,  DX  is  designated  as  an 
operator  for  differentiation.    The  scope  of  the  DX  operator,  i.e.  , 
its  variables  of  differentiation  and  its  expression  to  be 
differentiated,  are  assigned  by  enclosure  in  parentheses  in 
accordance  with  the  following  format:    DX;  the  left  parenthesis; 
the  list  of  variables  of  differentiation,  each  separated  by  a 
comma ,  and  the  list  terminated  by  a  double  period;  the 
expression  to  be  differentiated;  the  right  parenthesis.    For 
example,  the  algebraic  equation, 

n2  2    2  \  2 

A  =  X         O  (X  Y    +         O     (XY  SIN   Z)) , 

^  X    ^Y  ^  Z 

would  be  expressed  in  Fortran  63  as  follows: 

A  =  X  *  DX(X,Y.  .X**2  *  Y**2  +  DX(Z.  .X  *  Y  *  SIN(2)**2)) 


2.    Translation  of  the  Derivative  Dual-Operator  and  its 
argument  into  reverse  Polish  notation. 

The  stack  compilation  technique  described  by  Wegner      is 
utilized  in  the  transformation  process. 

In  order  to  determine  its  "scope"  in  a  reverse  Polish  string, 
the  DX  operator  must  assume  the  dual  nature  afforded  the 
Fortran  "Function".     Therefore,  the  DX  is  assigned  a  suffix 
operator;  namely,  a  double  comma  (the  obvious  "prime"  not  being 
available  on  a  keyboard) .    When  DX  is  encountered  during 
translation  from  Fortran  into  reverse  Polish,  control  passes  to  the 
DX  closed  subroutine  (Figure  4     Page  23)  which  performs  the 
following  functions.    The  "DX"  is  released  immediately  into  the 
Polish  string,  and  its  suffix  operator  is  placed  into  the  operator 
shunt  with  a  priority  assigned  sufficiently  high  to  insure  its 
release  immediately  upon  the  encounter  of  the  first  operator 
following  the  expression  being  differentiated.    The  argument's 
left  parenthesis  is  then  placed  into  the  shunt,  the  variables  of 
differentiation  extracted  and  placed  into  the  Polish  string,  the 


Wegner,  P.  ,  Introduction  to  System  Programming, 
Academic  Press,   1964:    100-121. 


double  period  placed  into  the  Polish  string,  and  control  then 
returned  to  the  superior  routine  for  normal  processing  of  the 
remainder  of  the  expression. 

3.     Execution  of  the  Derivative  Dual-Operator. 

Upon  encountering  the  DX  dual-operator  during  the 
translation  of  the  reverse  Polish  notation  into  a  machine 
symbolic  language,  control  passes  to  the  DX  closed  subroutine. 
(Figure  5) 

First,  the  variables  of  differentiation  are  extracted,  indexed, 
and  stored.    The  double  period  signifies  the  beginning  of  the 
expression  to  be  differentiated.    The  differentiation  process  is 
then  executed  with  respect  to  the  variables  in  the  order  listed  by 
the  programmer. 

Differentiation  is  executed  upon  operands,  whether  single 
symbols  or  expressions,  in  the  order  of  presentation  of  their 
algebraic  operators,  in  the  normal  left  to  right  scan.    A  simple 
algorithm    (Page  27  )  is  employed  to  identify  the  operands  of  each 
algebraic  operator  encountered.    The  operands  are  identified  with 
respect  to  position  and  complexity,  and  their  derivatives  are 
identified  as  zero  or  otherwise.    The  differentiation  process  is 
then  executed  in  accordance  with  the  applicable  rule  (Page  28  ). 
The  partial  result  thus  obtained  is  stored  in  an  array,  which 


serves  as  a  base  upon  which  to  build  the  derivative  of  the 
following  algebraic  operator  and  its  pair  of  operands.    In  most 
cases,  this  partial  result  becomes  the  first  part  of  the  derivative 
of  the  following  operation;  thus  symbol  manipulation  is 
minimized.    In  the  remaining  few  cases  only  part  of  the  partial 
result  requires  shifting.     Further,  zeroes  are  identified 
immediately,  eliminated,  and  the  process  scan  proceeds  on  to  the 
next  algebraic  operator. 

Upon  completion  of  execution  of  partial  differentiation,  the 
symbolic  result  is  patched  into  the  reverse  Polish  string  between 
the  DX  and  its  suffix.    At  this  point,  the  expression  differentiated 
and  its  derivative  can  be  printed  in  reverse  Polish  notation  or 
transformed  into  normal  form  prior  to  print-out.    Though  not 
provided  in  the  attached  model  program,  a  code  word  could 
easily  be  inserted  into  the  DX  format,  following  the  variables 
of  differentiation,  where  it  could  easily  be  identified  and  used 
to  command  a  print-out  of  the  symbolic  derivative  only  when  req- 
uired   by  the  programmer  for  analysis. 

If  a  DX  operator  is  encountered  during  differentiation,  the 
partially  processed  results  are  stored,  the  new  DX  operator  is 
executed,  the  stored  partial  results  retrieved,  and  the  process 
continued  from  the  point  of  interruption. 

8 


IV.     TEST  RESULTS 

Sample  runs  of  the  program  were  made  on  the  CDC  1604 
computer  at  the  U.  S.  Naval  Postgraduate  School.    Thirty  test 
equations  were  used,  the  last  of  which  being  the  equation 
presented  "by  Hanson. 

The  time  required  to  run  the  entire  program,  including  the 
printing  out  of  the  original  input  equation  in  Fortran,  the 
translation  of  the  Fortran  equation  to  reverse  Polish  notation, 
the  derivatives  of  the  required  portions  of  the  equation  in  reverse 
Polish  notation,  and  the  printing  out  of  the  macros  for  the 
evaluation  of  the  equation  was  five  minutes,  fifty-three  seconds. 
The  time  required  for  the  compilation  of  the  program  was  three 
minutes,  five  seconds. 

Additional  runs  of  the  program  were  made  to  determine  the 
time  for  the  internal  execution  of  the  program,  eliminating  all 
print-outs.    The  time  for  this  run  was  four  minutes.    The  time 
required  to  execute  the  program  giving  all  print-outs  with  the 
exception  of  the  mapros  was  four  minutes,  twenty-eight  seconds. 

The  sample  programming  runs  that  were  made  showed  that 
thirty  equations ,  some  of  which  required  the  derivative  of  a 

derivative,  required  two  minutes,  forty-eight  seconds. 

11 

Hanson,  J.  W.  ,  et  al.  ,  op.  cit.  ,   349-355. 
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APPENDIX 


LIST  OF  DIFFERENTIALS  PROVIDED  IN  MODEL  PROGRAM 


d  (a  *  u) 

=  a  *  du 

d  (u  +  ;v) 

=?  du  +  dv 

d  (u  *  v) 

=  u  *  dv  +  v  *  du 

d  (u  /  v) 

=  (v  *  du  -  u  *  dv)  /  v2 

d  (un) 

=  n  *  u11™1  *  du 

d(uv) 

=  v  *  uv=1  *  du  +  uv  *  ln(u)  *  dv 

d  (au) 

=  a11  *  In  (a)  *  du 

d  (uu) 

=uu  *  (1  +  ln(u))  *  du 

d  exp(u) 

=  exp(u)  *  du 

d  ln(u) 

=  du/u 

d  sin(u) 

=  cos(u)  *  du 

d  cos(U) 

=  -sin(u)  *  du 

d  tan(u) 

=  sec2  (u)  *  du 

d  sqrt(u), 

=  .5*du/  sqrt(u) 

Note  1.    Additional  functions  can  be  easily  added  to  routine 
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z  o  o  o 


cc  cc 

UJ  LU 


0£ 

G 


<  < 

£V  £V 

LU  LU 

a.  o_ 

o  o 


c 


LU 

CC 


LU 
> 


LU 
v. 
LU 
j_ 
U_ 


L> 

o 

u_ 

O 

O 


O 


3  X 

o  o 


LU  UJ  <   < 

_l  _J  I-  > 

O  L0   tO  — < 

z  z  z  cc 

— ■  •—  o  uj 

tO  to  U  Q 


z 
o 

—  LU  LU 

^  _I  Cu 

tO  CO  o 

LU  <  I 

CC  — '  _J 

d.  or:  < 

X  <  x 

UJ  >  O 


tO  tO 
UJ 

>  o 


c  —  ~ 

Z   I—  tO 

•—  <  to 

>  UJ 

_  QT 

or:  a 

UJ  x 

Q  UJ 


01 

'— 
to 

to 


c 

to  Q. 

cc 

C   LU 
t-  lO 


o  o 


a: 

LU 


< 
or 

LU  > 

Cl  LU 
O  CC 


UJ  LU 
CD   CD 


z  z 

err  or 

I—    !— 

tO  to 

LU  Z 

>  o 

I—  to 
<  tO 
>  LU 

—  or: 
cc  a 
lu  x 

Q  LU 
Lu   Lu 

o  o 

q  a 
z  z 

Lu1  UJ 


Z  <    O 
O  Q  ~ 


or 

C 

I— 
< 

or 

a 
o 

or 
o 
u_ 

v." 


< 

:-: 

z 

or 

h- 


UJ 

O 

or 

CO 


< 


tO  51 

to  < 

UJ  _J 

or 

a  u_ 

X  o 

UJ 

LU 

LU  > 

CO  -- 


[/) 


u.  or 

LU  UJ 

_J  o 
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tO  o 

LU  X 

or  or 

Ol 

X  Lu 

LU  O 

LU  LU 

Q  > 

tO  h- 
< 

1-  > 

X  -- 

to  or 

—.  lu 

or  o 


to 
a 
z 
< 

or 

UJ 

or  cl  to 

C  O  uj 

Lu        _J 

LU  CD 

f-  _J  < 

Z  LO  — 

x  z  or 
>-  < 
W  > 


z 
to 


X 

o 

u 
o 


Ll 

o 

-: 
o 


to  — • 
O  tO 

a  o 

—  a 

o 
z  o 

— •  z 
or  — < 

HOT  O 

tO  1—  z 

to  *— 

x  or 

to  x  l- 

— .  to  to 

_j  — < 

OJI 

a  o  to 

o_  — ■ 

uj  _i 

to  uj  c 

or  to  a 
uj  or 

>  LU  LU 

Ul    >  tO 

or  lu  or: 

or  lu 

o  > 

t—  s:  lu 

z  o  or 

—,  or 

u_  u_ 

z  o 

o  z 

—  o  z 
i-  —  c 
I— I  I—  _1 
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O  to  •— 
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Cl 
>- 
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(-   — i  Q 

Z  X  z 

LU   LU  LU 
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a 


Z  O  LU  _J  Z 

•-  z  z  z  z  z 


Z  LU          tO 
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Q  LU 

—  z       z 

z  z 
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*-<  — -   Q  LU 

Cl 

Q 

o 
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X 

> 

_J  KD  O  O  Q 
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c 

o 

z 

~ 

or  z  r—  — ■  or 

z  c  a  to  > 

a 

■— 

X   0_ 
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— 
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z 
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or 

o- 
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tO 

>   LU  UJ   LU   1— 

Figure  3 


13 


TRANSFORMATION'  OF  FORTRAN  63  ALGEBRAIC  EQUATION 
INTO  REVERSE  POLISH  NOTATION 


1A. 


JENTER  WITH  ALGEBRAIC  EQUATION 

STORED  IN  ALPHA, 

ONE  CHARACTER  PER  WORD, 

LEFT  JUSTIFIED, 

j 

\7 

INITIALIZE  COUNTERS 


ADVANCE  ALPHA  POSITION  CTR. 


YES 


v 


IS  ALPHA(I)  A  BLANK  SPACE  ? 


v 

820 


YES 
7^~ 


V 


NO 


IS  ALPHA (I)  AN  OPERATOR  ? 


v/ 


NO 


ADVANCE  ALPHA  POSITION  CTR,  I. 


ADVANCE  CHARACTR  CTR,  K. 
SET  CHARACTR (K)   =  ALPHA(I) 


\/ 


IS  THE  CHARACTR  CTR  GREATER 
THAN  8  ? 


YES 


ERROR  ROUTINE. 

SYMBOL  EXCEEDS  MAXIMUM  OF 

8 "ALLOWABLE  CHARACTERS. 


IS  ALPHA (I)  AN  OPERATOR  ? 


NO 


I ' ^ 

SET  CHARACTR  CTR,  K,   =1 

SET  CHARACTR(K)  =  ALPHA(I) 


i/ 


805 


EXIT 


Figure  4 
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V 


IS  ALPHA(I)  M  "*"  ? 


ZES 


ADVANCE  ALPHA  POSITION  CTR. 


NO 


>k. 


IS  ALPHA(I)  AN  "*«■  ? 


=- 


l£ 


NO 


REDUCE  ALPHA  POSITION  CTR. 


\v 


ROUTINE  FOR  THE  TRANSFORMATION) 
OF  THE  CURRENT  OPERATOR  , 
WHICH  IS  IDENTIFIED  BY  THE 
INDEX  "N"  . 


Figure  4  (Cont'd) 
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\z 


REDUCE  ALPHA  POSITION  CTR. 


V 


ADVANCE  POLISH  POSITION  CTR. 

ENCODE  CHARACTR(J),   J=T.K 
INTO  POLISH(L).  FORMAT  (8A1) 


IS  ALPHA(I)  A  BLANK  SPACE  ? 


± 


NO 


IS  ALPHA (I)  A  LEFT  PAREN  ? 


YES 


.  NO 


POLISH (L)  IS  A  SYMBOL, 
ADVANCE  PARITY  CTR. 


ARE  THERE  ANY  MORE  CHARACTERS  j  YES 

s 


IN  ALPHA  ? 


V 


NO 


ffiCORD  LE1\TGTH  -OF  POLISH  STRING 


i  825 


802 


\y 


EXIT 


Figure  4  (Cont'd) 
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... 


w 


:S  POLISH (L)    A  DUAL-OPERATOR?  j- 

— — — — a— — ^— — — ■ ■  ■  t  -  ■   -     >■->*—  n.w'i   .  r   i  ■  ■  i  .  -  ■  ■-' 


YES 


4 


ROUTINE  FOR  TRANSFORMATION 
OF  A  DUAL-OPERATOR. 


NO 


Te   ^TTOTT/T   *A    WT.T«mT^„            1  YES    -sj   ROUTINE  FOR  TRANSFORMATION 
IS  POLISH(L)   A  FUNCTION  ?       j ^  Qp  A  g^^ 


I  NO 


IS  POLISH (L)  A  SUBSCRIPTED  jj  YES  xj  ROUTINE  FOR  TRANSFORMATION 
VARIABLE  ?      (NOTE  1)  ^1  OF  A  SUBSCRIPTED  VARIABLE. 


\    NO 


^. — ~ — _.^_^_.- 


ERROR  ROUTINE* 
EITHER  THIS  SYMBOL  IS  AN 
UNDIMENSIONED  ARRAY,  OR  AN 
OPERATOR  IS  MISSING  AFTER 
THIS  SYMBOL. 


NOTE  1 .   FOR  SIMPLICITY,  SUBSCRIPTED  VARIABLES  ARE  NOT  CONSIDERED 
HEREIN.   THE  NECESSARY 'ROUTINE  SHOULD  BE  INSERTED  INTO 
PROGRAM  AT  THIS  POINT,.  GIVING  CONSIDERATION  TO  PARITY 
COUNTERS  AND  TESTS.. 


Figure  4  (Cont'd) 
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\/ 


CURRENT  OPERATOR  IS  A 
"+"  OR  A  "-" 

TEST  FOR  UNARY"  OPERATOR. 


v 


NO 


IS  PARITY  EQUAL  TO  ZERO-  ? 


WAS  THE  PREVIOUS  OPERATOR  AN   N0 


11 


YES 


"="'  OR  AN  "("  ? 


YES 


r 


IS  CURRENT'  OPERATOR  A  "-"  ? 


NO 


w 


YES 


CURRENT  OPERATOR  IS  A 
UNARY  "-"  .   (NOTE  1) 


PLACE  A  ZERO  IN  POLISH  STRING 
INCREASE  PARITY  CTR. 


850 


^@>-^| 


8or 


ERROR  ROUTINE. 
OPERAND  MISSING. 


"/ 


NOTE  t.  A  UNARY  "-"   IS  TREATED  AS   (  0  -  OPERAND  ). 


Figure  4  (Cont'd) 
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REDUCE  PARITY  CTR. 


\y 


IS  PARITY  CTR  EQUAL   TO  0  ? 


\S/  ™ 


j  IS  THE  PRIORITY  OF  THE  TOP 
I  OPERATOR  IN  THE  SHUNT  LESS 
1  THAN  THE  PRIORITY  OF  THE 
CURRENT  OPERATOR  ? 


TRANSFER  THE  TOP  OPERATOR 
I  FROM  THE  SHUNT  TO  THE 
POLISH  STRING. 


■— 


NO 


890 


YES 


en 


ERROR  ROUTINE. 
OPERAND  MISSING. 


\V 


|mcTj 


S3, 


PLACE  THE  CURRENT  OPERATOR 
INTO  THE  OPERATOR  SHUNT. 


Figure  4  (Cont'd) 


19 


CURRENT  OPERATOR  IS  A 
"*"  OR  A  "/" 

SET  PRIORITY  EQUAL  TO  5, 


CURRENT  OPERATOR  IS  AN  "(" 
SET  PRIORITY  EQUAL  TO  2. 


CURRENT  OPERATOR  IS  A  "**» 
SET  PRIORITY  EQUAL  TO  6. 


851 


vv 


CURRENT  OPERAND  IS  A  FUNCTION! 
SET'  PRIORITY  EQUAL  TO  6. 


. & 

jjINSi^RT  THE  FUNCTION  OPERATOR, 
"."  ,  INTO  THE  OPERATOR  SHUNT 

PLACE  FUNCTION'S  ARGUMENT'S 
LEFT  PAREN  INTO  THE  OPERATOR 
SHUNT  AND  SET  ITS  PRIORITY 
EQUAL  TO  2. 

ADVANCE  PAREN  PARITY  CTR. 


w 


^801J 


801 
■ire  4  (Cont'd) 


on 


CURRENT  OPERATOR  IS  AN  ")"  . 
REDUCE  PAREN  PARITY  CTR. 


/ 


IS  PAREN  PARITY  COUNTER 
LESS,  THAN  0  ? 


LEFT  PAREN  MISSING 


il 


> 


V 


NO 


" 


YES  >vl  ERROR  ROUTINE.     t— ^JeXET 


IS  THE  TOP  OPERATOR  IN  SHUNT  '  YES  ^» 

A  LEFT  PAREN  ?       f— >>j}  REMOVE  LEFT  PAREN  FROM  SHUNT. 


NO 


TRANSFER  TOP  OPERATOR  FROM 
SHUNT  TO  POLISH  STRING. 


801 


/ 


i 


CURRENT  OPERATOR  IS  AN  "="  . 
SET  PRIORITY  EQUAL  TO  1 . 
REDUCE  PARITY" CTR. 


V 


Lis  polish  string  ctr,  l, 

EQUAL  TO  1  ? 


\7 


YES 


a 


"IS  SHUNT  CTR,  M,  EQUAL  TO  0  ?! 


NO 


3  853 


YES 


ERROR  ROUTINE. 
|    EQUAL   SIGN  OUT  OF  ORDER, 


\v 


?EXIT 


Figure  4  (Cont'd) 
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\!/ 


CURRENT  OPERATOR  IS  A  "$n 


V/ 


REDUCE  PARITY  CTR 


\y 


IS  PARITY  CTR  EQUAL  TO  0  ? 


W 


YES 


w 


TRANSFER  THE  TO?  OPERATOR 
FROM  THE  OPERATOR  SHUNT  TO 
THE  REVERSE  POLISH  STRING. 


NO 


IS  PARSN  PARITY  CTR 

EQUAL  TO  0  ?' 

t.  .  

,  v/ 

*\  IS  OPERATOR  SHUNT  EMPTY"  ? 
I . 


NO 


YES 


ERROR  ROUTINE. 

OPERAND  OR  OPERATOR  MISSING. 


h 


ERROR  ROUTINE. 
UNMATCHED  PARENS. 


RECORD  THE  LENGTH  OF  THE 
POLISH  STRING. 


V 


KXIT 


\7 


Figure  4  (Cont'd) 
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POLISE(L)  IS  EQUAL  TO  THE 
DERIVATIVE  DUAL -OPERATOR,  DX. 


V/ 


ADVANCE  OPERATOR  SHUNT  CTR,  M. 


INSERT  THE  DUAL-OPERATOR 
SUFFIX  ",,"  INTO  SHUNT. 

SET  PRIORITY  EQUAL  TO  6. 


w 


ADVANCE  OPERATOR  SHUNT  CTR,  M, 

PLACE  ARGUMENTS  LEFT  PAREN 
INTO  THE  OPERATOR  SHUNT,  WITH 
PRIORITY  EQUAL  TO  2. 

INCREASE  PAREN  PARITY" CTR. 


V 


ADVANCE  ALPHA  POSITION  CTR,  I.  <$ 


NOTE  1 .  FOR  SIMPLICITY,  THIS  ROUTINE  RESTRICTS  VARIABLES  OF 
DIFFERENTIATION  TO  A  SINGLE  CHARACTER  SYMBOL.   IF 
DESIRED,  THE  ROUTINE  PREVIOUSLY' DESCRIBED  FOR 
.  ASSEMBLY  OF  A  MULTIPLE  CHARACTER  SYMBOL  SHOULD  BE 
INSERTED  AT  JUNCTION  882. 


Figure  4  (Cont'd) 
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t 3 \ 

\   ADVANCE  POLISH  POSITION  CTR,  L 

SET  POLISH(L)  EQUAL  TO  ALPHA(I) 


V 


ADVANCE  ALPHA  POSITION  CTR,  I  l<$ 


V 


IS  ALPHA(I)  A  BLANK  SPACE  ? 


YES 


W 


NO 


IS  ALPHA(I)  A  COMMA  ? 


'■"' 


NO 


IS  'aLPEA(I)  A  PERIOD  ? 


NO 
I— $» 


^ 


YES 


ADVANCE  ALPHA  POSITION  CTR,  I. 
IS  ALPHA(I)  A  PERIOD? 


NO 


■■'/ 


YES 


ADVANCE  POLISH  POSITION  CTR5  LJ 
SET  POLISH (L)  EQUAL  TO.  ".." 


801 


ERROR  ROUTINE. 

DOUBLE  PERIOD  MISSING. 


^ :  ^ 


EXCT: 


Figure  4  (Cont'd) 
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DERIVATIVE  ROUTINE 


© 

v/ 


rn.jT  v 


THE  CURRENT  POLISH  CHARACTER 
IS  THE  DUAL-OPERATOR  DX. 


vv 


SET'  THE  DERIVATIVE  CTR,  ND? 
EQUAL  TO  ZERO* 


0 


IS  THE  DERIVATIVE  CTR 
GREATER  THAN  ZERO? 


YES 


> 


TRANSFER  CURRENT  DATA 
TO  SCRATCH. 


v 


NO 


SET  POLISH  POSITION  CTR,  NP, 
.  EQUAL  TO  "ENTRANCE"- 


w 


ADVANCE  THE  DERIVATIVE  CTR.. 


<c 


v/ 


SET'  CTRS  ND,  NE,  NL,  NS,  NV 
EQUAL  TO  ZERO;  AND 
NX  EQUAL  TO  ONE. 


\v 


SET  "ENTRANCE"  EQUAL  TO 
POLISH  POSITION  CTR. 


'ADVANCE  POLISH  POSITION  CTR. 


VV 


ADVANCE  VARIABLES  CTR. 
JJSET  VARIABLE (NV)  EQUAL  TO 
| POLISH  (NP). 


\/ 


; 


'  IS  POLISH(NP)   A  DOUBLE  PERIOD? 


NO 


Figure  5 


\  / 


k^3 
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^ 


SET  THE  VARIABLES  INDEX,  NIV, j 
EQUAL  TO  NX. 


/ 


6-1- 


\/ 


SET  THE  CURRENT  VARIABLE,  V, 
EQUAL  TO  VARS(NIV).. 


|SET  THE  OPERANDS  CTR,  NO, 
EQUAL  TO  ZERO. 


> 


ADVANCE  POLISH  POSITION 
CTR,  NP. 


^S  -IS  POLISH(NP)  A  ELANK  SPACE  ? 


± 


\7 


NO 


|  IS  POLISH(NP)  AN  OPERATOR  ? 


YES 


v 


1/  NO 


IS  POLISE(NP)  A  DUAL-OPERATOR?! 


YES 


11 


_2_ 
■^1600 


\y 


ADVANCE  THE  OPERANDS  CTR,  NO. 


\7 


.ADVANCE  OPERAND  SHUNT  CTR,    NS, 
ji  PLACE  POLISE(NP)    INTO   TEE 
OPERAND  SHUNT. 


Figure  5  (Cont'd) 
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2  OR  MORE 


£00) 


© 

\iy 

CLASSIFY  TEE  NUMBER  OF 
OPERANDS  IN  THE  CURRENT 
UNBROKEN  SEQUENCE. 

j  DIE 
% 


\/ 


-- 


j ADVANCE  SEQUENCE  COUNTER,  NL, 


\/ 


I  THE  CURRENT  OPERATOR  IS 
H APPLICABLE  TO  THE  LAST  TWO 
I  OPERANDS  IN  CURRENT  SEQUENCE. 
i 

i REDUCE  THE  COUNT  OF  THE 
j OPERANDS  IN  THIS  SEQUENCE  BY' 
2,  AND  RECORD  THE  COUNT.. 


W 


PROCEED  TO  ROUTINE  FOR 
DIFFERENTIATION  OF  TWO 
SINGLE  OPERATORS. 


tT5T) 


~_ 


ZERO 


(300) 


HAVE  ALL  OPERANDS  IN  THE  LAST! 
SEQUENCE  BEEN  PROCESSED  ? 


YE 


w 


NO 


THE  CURRENT  OPERATOR  IS 
APPLICABLE  TO  THE  LAST 
SEQUENCE. 

REDUCE  THE  COUNT  OF 
UNPROCESSED  OPERANDS  IN 
THAT  SEQUENCE  BY  1.. 


\7 


PROCEED  TO  ROUTINE  FOR 
DIFFERENTIATION  OF  A 
LEFT  SIDE  SINGLE  OPERAND  AND 
A  RIGHT  SIDE  EXPRESSION. 


v/  m 


1  M.P>iHi. 


\V 


i CURRENT  OPERATOR  IS  APPLIC- 

THE  CURRENT  SINGLE 
AND  THE  PROCESSED 
iPORTlON  OF  THE  PRECEEDING 


ABLE  TO 
^OPERAND 


■SEQUENCE. 


y 


IpOCEED  TO  ROUTINE  FOR 
[DIFFERENTIATION  OF  A  LEFT 
[SIDE  EXPRESSION  AND  A  RIGHT 
(SIDE  SINGLE  OPERAND,, 


vy 


200 


C      /r,nni-tA\ 


TL 


CUKRENT  OPERATOR  IS  APPLIC- 
ABLE TO  THE  LAST  SEQUENCE  AND 
THE  PROCESSED  PORTION  OF  THE 
PRECEEDING  SEQUENCE. 


v/ 


REDUCE  THE  SEQUENCE  CTR,  NL. 

— H 


PROCEED  TO  ROUTINE  FOR 
DIFFERENTIATION  OF  WO 
EXPRESSIONS. 


v 


■ 

CM 

V 

1 

O 

~> 

e> 

s 

IH 

CD 

Figure  5  (Cont'd) 
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Figure  5  (Cont'( 


2$ 


<*■> 

03 

S 

^-> 

H 

D 

1 

k          - 

Eg 

CO 

1     s 
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Figure  5  (Contad) 
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CO 

a 
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Figure  5  (Cont'd) 
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^i    VI 


. 


Figure  5  (Cont'd) 
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w 


POLISH(KP)  IS  THE 
DUAL-OPERATOR  SUFFIX  ",,", 


PRINT  THE  EXPRESSION  AND  ITS 
SYMBOLIC  DERIVATIVE  IN  REVERSE 
POLISH  NOTATION*   (OR  RETURN- 
JUMP  TO  ROUTINE  TO  CONVERT  TO 
NORMAL  FORTRAN  NOTATION*) 


AZ 


'SET  "EXITU  EQUAL  TO  THE 
OLISH  POSITION  COUNTER. 


.TCH  THE  DERIVATIVE  INTO  THE 
REVERSE  POLISH  STRING  BETWEEN 
|THE  "DX"  AND  THE  ",,"* 


\/ 


j. 


|SET  THE  POLISH  POSITION  CTR 
EQUAL  TO  "ENTRANCE". 

SET  CTRS  ND,  NE,  NL,  NS 
EQUAL  TO  ZERO. 

ADVANCE  VARIABLES  INDEX,  NIV. 


(HAS  DIFFERENTIATION  BEEN 
COMPLETED  FOR  EACH  VARIABLE? 


Sno 


w 


YES 


PSLETE  THE  DUAL-OPERATOR  FROM 
THE  REVERSE  POLISH  STRING. 


■c**« 


> REDUCE  THE  DERIVATIVE  CTR. 


NO 


S  THE  DERIVATIVE  ARRAY  EMPTY? — » 

r 


\i/_ 


YES 


1 

THE  EXPRESSION  BEING  DIFFERENT 
TIATED  CONSISTS  OF  A  SINGLE 
SYMBOL*   TEST1  THE  SYMBOL  AND 
THEN  SET  THE  DERIVATIVE  EQUAL 
TO  "ONE"  OR  "ZERO". 


ERROR  ROUTINE, 


IS  THERE  AN  INTERRUPTED 
"j  DERIVATIVE  OPERATOR,  DX*.  ? 


|   YES 


BACKSPACE  THE  SCRATCH  TAPE. 

TRANSFER  DATE  FROM  SCRATCH. 

BACKSPACE  THE  SCRATCH  TAPE. 

SET"  NX  EQUAL  TO  VARIABLES 
INDEX. 


-  — 

6v 


Figure  5  (Cont'd) 
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